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Abstract 

We demonstrate by simple mathematical considerations that a power-law tailed distribution in the 
kinetic energy of relativistic particles can be a limiting distribution seen in relativistic heavy ion experi- 
ments. We prove that the infinite repetition of an arbitrary composition rule on an infinitesimal amount 
leads to a rule with a formal logarithm. As a consequence the stationary distribution of energy in the 
thermodynamical limit follows the composed function of the Boltzmann-Gibbs exponential with this for- 
mal logarithm. In particular, interactions described as solely functions of the relative four-momentum 
squared lead to kinetic energy distributions of the Tsallis-Pareto (cut power-law) form in the high energy 
limit. 

1 Why to investigate composition rules 

One of the theoretically challenging questions related to relativistic heavy ion physics is how to 
establish the existence of a thermal or near thermal state of (quark) matter during high energy 
collisions. Besides studies of the hadronic flavor composition (TJ [2 [H [5] , the shape and the 
steepness of transverse momentum spectra on different particles presents experimental informa- 
tion on this question [SI [7J [HI M EH EJ HI]. Since f° r high energies of the out-coming particles 
(pions, kaons, antiprotons, etc.) these spectra show a power-law tail, it is important to clarify, 
whether this prominent feature can be explained in the framework of general thermodynamic 
ideas, or it has to be rendered to the realm of case by case non-equilibrium dynamics. While in 
some earlier works we have demonstrated that such spectra can be obtained as stationary distri- 
butions by altering the energy composition rule in two-particle collisions from the simple addition 
to another rule [HI [13] , in the present paper we aim at understanding the general mechanism 
setting an effective rule and a non-exponential stationary distribution of the individual particle 
energies in the thermodynamical limit. 

The classical thermodynamics using extensive and intensive quantities is being extremely suc- 
cessful in describing and understanding uncountable physical phenomena in nature. However, 
there are particular cases, where the observed distribution of individual energies in a complex 
system does not follow the exponential Boltzmann-Gibbs law. There are suggestions trying to go 
beyond the classical picture and to consider descriptions generalizing traditional thermodynam- 
ical concepts, in particular to consider the possibility that the composition rule deviates from 
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the simple addition. Such quantities are often called 'non-extensive', although strictly speaking 
the extensivity property is required only in the thermodynamical limit, i.e. for large systems 
consisting of many particles. 

In particular generalizations of the entropy formula, connecting the quantity of macroscopic 
entropy to probabilities of microstates of an extended system, have been repeatedly suggested 
in forms generalizing the Boltzmann-Gibbs-Shannon logarithmic formula [T51 [TBI [T71 [T5] . The 
deformed (or generalized, extended) logarithm relates the abstract product rule to a summation 
or composition formula in a general way: the statistical independence of states is hence mapped 
to a non-additivity of the composition formula[3U [551 H31 [Ml US]. Reversely, the additivity of 
entropy is - in some cases - achieved by non-product probabilities, trying to grasp the essence of 
surviving correlations (surmised to occur due to long range interactions) in systems, which are 
large in the thermodynamical sense [26]. The additivity of entropy also can be achieved by weakly 
non-local extensions [19l [20] . Naturally, the use of a formal composition rule is also applicable to 
the energy [27l [28]. Two subsystems combine to a common larger system not necessarily additively 
with their energies; the interaction part may lead to a finite relative contribution in the infinite 
particle number limit. 

2 Repeated compositions and the proof of associativity 

The important question arises, that what happens if we repeat (compose with itself) an abstract 
mathematical composition rule (the formal product in group theory) infinitely many times, but 
each time applying to an infinitesimal amount: grasping this way the extensivity limit of arbi- 
trary composition rules. We find that some rules transform to the simple addition in this limit, 
while others not. The proposition is that the thermodynamical limit of an arbitrary pairwise, 
iterable composition rule is an associative rule. The importance of this statement becomes clear 
by referring to the known mathematical property that associative rules always possess a strict 
monotonic function, called here the formal logarithm, in terms of which they can be expressed |29j. 
Let us denote an abstract pairwise composition rule by the mapping (x,y) — > h(x,y). Whenever 
h(x,y) is an element of the same set as x and y, the composition is iterable arbitrarily long. 
Associativity of such a rule is formulated by the function equation 

h(h(x,y),z) = h(x,h(y,z)) (1) 

for x, y and z being elements of the same group. For our purpose we shall consider energies or 
entropies of physical subsystems. The general solution of the associativity equation (UJ is given 

by 

h(x,y)=X- 1 (X(x)+X(y)) (2) 

with X(x) being a strict monotonic function. It is referred to as the "formal logarithm" , because 
it maps the arbitrary composition rule h(x, y) to the addition by taking the X-function of eq.©: 

X(h(x,y))=X(x) + X(y). (3) 

Due to this construction the generalized analogs to classical extensive (and additive) quantities 
are formal logarithms, whenever the composition rule is associative. As a consequence station- 
ary distributions, in particular by solving generalized Boltzmann equations [TJ], are the Gibbs 
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exponentials of the formal logarithm, 



/(*) = ie-W>. ( 4 ) 



Now we investigate a large number of iterations, N, of the composition rule applied to an 
infinitesimal amount y/Nm each step: 

x n(v) ~^°-y °h(j^>---,jf) ■ (5) 
N-l 

Whenever the limit, 

lim xn(v) < oo, (6) 

N— »oo 

is finite for a finite y, we are dealing with an extensive (but not necessarily additive) system. Our 
purpose is to study such systems and to obtain their asymptotic composition rule, 

XNt+Ns = <p(%Ni,XN 3 ) (7) 

in the limit N±,N 2 — > 00. The repetitive composition can be formulated as a recursion at an 
arbitrary step n between and N as follows: 

x n = h (x n -i, jj) , (8) 

with xq = 0. It is a natural requirement, but important for that what follows, to consider only 
such rules which satisfy h(x, 0) = x. Subtracting x n -i = h(x n -x, 0) from this formula we arrive 
at 

y 

x n - x n -i = h(x n -i, — ) - h(x n -i,0). (9) 

Introducing the (in the thermodynamical limit continuous) variable t — ne, we follow an evolution 
in t alike the renormalization flow: 

In the expression on the right hand side h' 2 (x, + ) denotes the partial derivative of the rule h(x, y) 
with respect to its second argument taken at this argument value approaching zero from above. 
The final time is given by tj = Ne. Note that the uniformity of subdivisions to y/N is not really 
necessary; all infinitesimal divisions summing up to tf lead to the same differential flow equation. 
The solution of eq. (fT0|) is given by 

X 

L{x)= / = y-. 11 

J h' 2 {z,0+) t f 


This solution, when strict monotonic and hence invertible, defines the following asymptotic com- 
position rule in the thermodynamic limit: 

x 12 := (f(x 1 ,x 2 ) = L^ 1 (L(ari) + L(x 2 )) , (12) 

due to L{x\) = yti/tf, L(x 2 ) — yh/tf and L(xi 2 ) — y{t\ +fa)/tf- Note that U/tf — Ni/N 
are the extensivity shares of the respective subsystems. This asymptotic composition rule is 
associative and commutative. 
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3 Classification, important examples 

Now we turn to the analysis of important particular rules and their asymptotic pendants in the 
thermodynamic limit. 

The trivial (and classical) addition is the simplest composition rule: h(x, y) = x + y. In this 
case h' 2 (x, + ) = 1 and one obtains 



L(x) = J dz = x, (13) 
o 

and with that the original Gibbs exponentials, e~@ E /Z, for stationary distributions of any Monte 
Carlo type algorithm using the original composition rule. The corresponding asymptotic rule is 
(p(x,y) =x + y. 

The rule leading to the g-exponential:2h (or Pareto, or Tsallis) distribution is given by 
h(x,y) = x + y + axy with the parameter a proportional to q — 1. In this case one obtains 
h' 2 (x, + ) — 1 + ax and 

X 

L(x) = / — — — = - lnfl + ax). (14) 
J 1 + az a 
o 

This formal logarithm leads to a stationary distribution with power-law tail as the function 
composition exp o L on the power — /3: 

f(E) = |e-i = \ (1 + aE)-^ a . (15) 

A generalized entropy formula on the other hand can be constructed as the expectation value of 
the inverse of this function, L^ 1 o In: 

S= f = -/(/ 1 - Q -/)- (16) 



The asymptotic composition rule again coincides with the original one: <p(x, y) — x + y + axy. 

A simple rule suggested by Kaniadakis[3T] is based on the sinh function. The formal logarithm 
is given as 

L(x) = -Arshfra), (17) 

and its inverse becomes i _1 (t) = smh(nt) / k. The stationary distribution, composed by exp o L, 
is 

/eq(?) = I («p + Vi + «V) _/3/K ■ (18) 

For large arguments it gives a power-law in the momentum absolute value p and hence also in 
the relativistic energy. The corresponding entropy formula is the average of L^ 1 o In over the 
allowed phase space: 

isinh(Kln/)= / J - ^ . (19) 
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The composition formula can be reduced to 



h(x, y) = x\J\ + K 2 y 2 +yy/l + k 2 x 2 . (20) 

For low arguments it is additive, h(x, y) s=s x + y, for high ones it is multiplicative, h(x, y) k, 2nxy. 
It has been motivated by the relativistic kinematics of massive particles. Regarding k = \/mc, 
up = sinh?7, so the formal logarithm is proportional to the rapidity, L(p) = mcrj. This implies a 
stationary distribution like exp{— 0mc\r]\), which has not yet been observed in relativistic particle 
systems. Therefore it is wishful to consider some further scenarios based on other quantities 
deduced from relativistic kinematics (see next section). 

The rule leading to a stretched exponential stationary distribution is given by h(x, y) = 



y b ) 1/b . Here some care has to be taken, the partial derivative has to be evaluated not at 
zero, but at a small positive argument, e = y/2N. We get h' 2 (x,e) = c(e)x 1 ~ b with a factor 
depending on e and - depending on b - possibly diverging in the e = limit. However, this 
does not spoil our procedure; we obtain L(x) = c(e)x b /b and with that the asymptotic rule: 

ip(x,y) — (x b + y 6 ) • The reason is that constant factors in the formal logarithm can be 
eliminated without loss of any information. 

Our next example is a non-associative rule; its asymptotic pendant cannot be itself. We 
consider 

y) = x + y + a — ■ — (21) 
x + y 

(a combination of arithmetic and harmonic means). The fiducial derivative is given by h' 2 {x, + ) = 
1 + a and - being a constant - it leads to L(x) = x/(l + a) and with that to the addition as 
asymptotic rule: ip(x, y) — x + y. 

As a last example in this train we discuss Einstein's formula for velocity addition, 

^ = TT^' (22) 

This rule is associative, and it also preserves its form in the thermodynamic limit. The fiducial 
derivative is given by h' 2 (x, + ) = 1 — x 2 /c 2 and the formal logarithm, L(x) = catanh(x/c) turns 
out to be the rapidity. The asymptotic composition rule recovers the original one. 

In the case of a general second order polynomial for the fiducial derivative h' 2 (z, 0) the asymp- 
totic composition rule turns out to be: 

x + y + axy 

f[x,y) = —— (23) 

1 + xy/c z 

with c 2 = — ziz 2 and a = — {z\ + z 2 )jz\z 2 \ Z\ and z 2 being the algebraic roots of h' 2 {z,0). It is 
a generalization of the Tsallis and Einstein rules. A similar composition rule has been found for 
the parallel transmittivity in certain Potts models, for a review see Rcf. 32 . 

Finally we prove that all associative rules are mapped to themselves by the thermodynamic 
limit. Given an original composition rule, h(x, y), which is associative, it can be expressed by its 
formal logarithm: 

h(x,y)=X- 1 (X(x)+X(y)). (24) 
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Then letting to act the strict monotonic function X on both sides and derive with respect to the 
second argument we obtain X'(h) dh/dy = X'(y) and 



h' 2 (x,0+)= *'® y (25) 
X'(h{x,0)) 

Due to the property h(x, 0) = x (equivalently X(0) = 0) the formal logarithm of the asymptotic 
composition rule is given by 

L { x)=[^d Z =^; (26) 

1 ' J x'(o) x'(oy K ' 

o 

it is proportional to the formal logarithm of the starting rule. Therefore the asymptotic rule is 
the same as we begun with: ip(x,y) — h(x,y). Actually the freedom in a factor of the formal 
logarithm always can be used to set X'(0) = 1. 

This way any associative composition rule describes a thermodynamical limit of a class of 
non-associative rules. Associativity is synonym to the thermodynamical limit. 

In Fig[T] we show the composition of y = 1 from y/2iV-sized pieces by the rule h(x,y) = 
x + y + G(xy) for G(w) = aw/(l + bw) with different parameters a and 6. While the 6 = 
case represents the Tsallis-Pareto rule (top figure, a), which asymptotically establishes, the 6 = 5 
case plot deviations from this for finite N-s (bottom figure, b). The 21-st point is the asymptotic 
composition of half-sized systems, </?(xjv , x/v), the other points plot x n by the recursive application 
of the rule. One inspects some deviation for 6^0, while for an associative rule h{x, y) no deviation 
in the end result occurs. 

Fig[2] shows the approach to the asymptotic limit as a function of the number of repetitions, 
N of the micro-rule, h(x,y) on a logarithmic scale. The series of X2N values is constructed 
by using h{x,y) = x + y + xy/(l + 5xy) - which is not an associative rule. The asymptotic 
rule corresponding to this choice is given by ip(x,y) = x + y + xy. The higher points belong 
to ip(xn,xn), they deviate from X2N f° r finite N. The continuous curve represents the 2iV-fold 
composition of 1/2N according to the asymptotic rule tp(x,y). It is given by the formula 

n / i ^ 2N 




X2N = < po ...c l p [—,...,— ) = + — ) -l. (27) 

The asymptotic rule is occurring slowly in this case, the convergence rate is given by the conver- 
gence of the Euler formula to the Euler number (xoa = e — 1). 

4 Deriving composition rules from interaction energy and 
kinematics 

In this section we review a few general considerations which may relate the interaction energy 
due to some kind of correlation to the use of formal logarithms and abstract composition rules. 
Our basic assumption is that the interaction energy between two subsystems ('particles') can 
be expressed as a function of the individual energies without the interaction (in asymptotic free 
states). This way E12 = E\ + E2 + U{E\, E2) is a general energy composition rule. 
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Numerical repetation of h(x,y) = x + y + axy 
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Figure 1: The series of composed amounts xn by subdividing y = 1 to N parts. In the top figure (a) 
results with the asymptotic rule h{x,y) = x + y + axy, while in the bottom (b) with a non-associative rule 
h(x,y) = x + y + xy/(l + 5xy). The 21-st values are the asymptotic compositions of half systems: y(ario,#io)- 



Convergence of h(x,y) = x + y + xy/(1 + 5xy) to the asymptotic rule 
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Figure 2: Approach to the asymptotic extensivity rule as a function of the number of repetitions of the 
non-associative rule, h(x,y) — x + y + xy/(l + 5xy), to the amount 1/N (boxes). Plotted are X2N arid for 
comparison the amount <p(xn , xn) as well as the curve corresponding to the 2A^-fold composition of if. 
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In most cases, discussed in physics, the interaction is given as a function of the relative 
distance, and this form is not directly related to kinematic data. In the sense of medium long 
time behavior, however, the interaction energy can often be expressed by the relative momentum: 
either due to a virial theorem or due to direct quantum mechanical solution for the relative wave 
function, like e.g. in the two-body Coulomb problem. 

In the followings we shall assume that the interaction energy is a function of the kinematic 
variable Q 2 , the square of relative four- momentum. We shall study whether relativistic speeds 
alone can cause " non-extensivity" , i.e. a power-law tailed kinetic energy distribution. The 
relativistic formula for Q 2 is given by the following Lorentz-invariant quantity: 

Q 2 = (pi - vi? ~ (Ei - E 2 f (28) 

with Pi,Ei being relativistic momenta and full energies of interacting bodies. Expressed by the 
energies and the angle O between the two momenta this becomes a linear expression of cosG: 

Q 2 = 2 (EiE 2 -m 2 -P1P2 cos 9) (29) 



with pi = E 2 — m 2 for i — 1,2. Here we use relativistic units (c = 1) and assume the same 
mass for both interacting partners, for simplicity. It is useful to note that writing eq. (129|) as 
Q 2 = 2(A - B cos 6) we have 

A ± B = E X E 2 - m 2 ± p lP2 = m 2 (cosh^i ± 772) - 1) (30) 

with using the rapidities to. For equal momenta or rapidities Q 2 = and A = B; A 2 — B 2 = 
m 2 (Ei — E2) 2 measures the energy difference. 

In order to estimate the interaction contribution we subtract the zero momentum terms, and 
assume 

E 12 = Et + E 2 + U(Q 2 ) - U(Ql) - U{Ql) + 17(0), (31) 

with Q 2 = 2m(Ei — to). This construction ensures that for equal momenta, i.e. for Q 2 — 0, no 
interaction correction occurs to the addition law for the energy. 

Seeking for an effective energy composition rule as an isotropic average over the relative 
directions of the respective momenta, one averages over the angle 9: 



(U(Q 2 )) = i / U(2A-2B cose) sinQ dQ 


F(2A + 2B) - F{2A - 2B) 



AB 



(32) 



with U (w) — dF/dw. This is easy to see upon the substitution w = 2(A — B cos 9). Since at zero 
momenta the total relativistic energies are nonzero, it is more physical to consider the kinetic 
energy only, Ki — Ei — m. The rule for the kinetic energy composition is hence given by 

F(2A + 2B)-F(2A~2B) 

Al2 = A1+A2H J75 

-U(2mK x ) - U(2mK 2 ) + 17(0). (33) 
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The coefficients A and B are also expressed by the respective kinetic energies: 



A = m{K l +K 2 ) + K l K 2 , 

B = K 1 K 2 {l + 2m/K 1 ) 1/2 {l + 2m/K 2 ) 1/2 . (34) 

One observes that the product of kinetic energies occurs due to kinematic reasons. 

Taylor expanding the integral of the unknown function U (w) around w = 2 A we obtain the 
following composition rule for the relativistic kinetic energies: 

h(x,y) = x + y - U{2mx) - U(2my) + U(0) 

+ E/ 2 ^)|^y (35) 

with A = m(x + y) + xy and AB 2 = 4xy(x + 2m)(y + 2m). The fiducial derivative becomes an 
expression with a finite number of terms 



h' 2 {x,0) = l-2mU'(0) + 2(m + x)U'(2mx) 
4 

+ -mx(2m + x)U"(2mx). (36) 



For all traditional approaches the interaction energy U is not considered as dependent on Q 2 . 
In such cases h! 2 {x,Q) = 1 and one arrives at the simple addition as composition rule. As a 
consequence the stationary energy distribution is of Boltzmann-Gibbs type. For Q 2 dependent 
interaction it is enlightening to analyze two particular kinematical cases: the extreme relativistic 
and the non-relativistic ones. In the first case m = has to be replaced and one obtains 

ti 2 (x,0) = l + 2xZ7'(0). (37) 

As discussed above this leads to a Tsallis-Pareto distribution in the kinetic energy (at zero rapidity 
in the variable % ~ m )- The opposite extreme, m 3> x leads to an undetermined asymptotic 
composition rule due to 

h' 2 (x, 0) w 1 + 2m (U'(2mx) - U'{0)) + ^m 2 xU"{2mx). (38) 

This result includes for U' = the traditional momentum independent interaction case leading to 
the addition as asymptotic rule for non-relativistic kinetic energies, and hence to the Boltzmann- 
Gibbs distribution. We note that in the relativistic kinematics the next simplest assumption, 
U' = a = const, leads to a Tsallis-Pareto distribution due to h' 2 (x, 0) = 1 + 2ax form eq . ([36f . 
while its non-relativistic approximation form eq. (|38p is still the Boltzmann-Gibbs exponential. 



5 Conclusion 

Summarizing we have proved that in the thermodynamical limit, by composing a finite amount 
of an extensive physical quantity as a repeated composition of infinitesimally small amounts of 
the same quantity one arrives at an effective asymptotic composition rule which is derivable from 
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a formal logarithm. This formal logarithm, L is constructed from the original composition rule 
uniquely and serves as a basis of describing a stationary distribution, expoL and the corresponding 
formal expression for the entropy it canonically maximizes, (L _1 o In). Associative rules lead to 
themselves in this limit, and are attractors for other rules. 

The addition is the simplest composition rule, the formal logarithm being the identity map. 
It leads to the Boltzmann-Gibbs distribution. The next simplest one leads to the Pareto-Tsallis 
power-law tailed stationary distribution. We also have shown that considering interaction ener- 
gies dependent on the relative four-momentum, in general a nontrivial asymptotic rule arises. In 
particular the Pareto-Tsallis distribution - and the corresponding non-additive asymptotic com- 
position rule - emerges generally from extreme rclativistic kinematics. In fact high energy particle 
spectra often show a power-law like tail in their kinetic energy. Besides that also non-relativistic 
systems may show such thermodynamical limit, but here the general case can also be different. 
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